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We consider the phenomenon of mutual synchronization in a fundamental quantum system of two 
detuned quantum harmonic oscillators dissipating into the environment. We identify the conditions 
leading to this spontaneous phenomenon showing that the ability of the system to synchronize is 
related to the existence of disparate decay rates and is accompanied by robust quantum discord and 
mutual information between the oscillators, preventing the leak of information from the system. 
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I. INTRODUCTION 



Synchronization phenomena have been observed in a 
broad range of physical, chemical and biological systems 
under a variety of circumstances In some instances, 
synchronization (also known as entrainment in this con- 
text) is induced by the presence of an external forcing 
or driving that acts as a pacemaker; in others it appears 
spontaneously as a consequence of the interaction be- 
tween the elements. The latter case is the most relevant 
from the point of view of complexity since it appears as 
an emergent phenomenon that takes place despite the 
natural differences between the elements. Collective syn- 
chronization, whose simplest description can be given in 
terms of coupled self-sustained oscillators, is found in re- 
laxation oscillator circuits, networks of neurons, cardiac 
pacemaker cells or fireflies that flash in unison Q. A 
key ingredient for collective synchronization is dissipa- 
tion, which is responsible for collapsing any trajectory of 
the system in phase space in a lower dimensional mani- 
fold. 



Synchronization has also been studied in the quantum 
world in the case of entrainment induced by an exter- 
nal driving Difficulties in addressing quantum col- 
lective synchronization come from the fact that in linear 
oscillators dissipation will lead to the death of the oscilla- 
tions after a transient while the extension to the quantum 
world of nonlinear phase oscillators does not allow for 
an insightful treatment. In this work, however, we take 
a step toward the understanding of quantum collective 
synchronization, showing that it is possible to fully char- 
acterize synchronization during the transient dynamics. 
We find that synchronization happens in the presence 
of a common bath due to a separation between dissi- 
pation rates. Different groups have recently approached 
this subject considering the (classical) synchronization of 
nanoscopic and microscopic systems susceptible of hav- 
ing quantum behavior [J-[6j . In this work we establish the 
connections between the phenomenon of synchronization 
and quantum correlations in the system. 



II. SYSTEM 

We consider two coupled quantum harmonic oscilla- 
tors dissipating into the environment with different 
frequencies jlQj | ? which is arguably one of the most fun- 
damental prototypical models. Current experimental re- 
alizations in the quantum regime include nanoelectrome- 
chanical structures and optomechanical devices pjj as 
well as separately trapped ions whose direct coupling has 
been reported recently [l2| . The system Hamiltonian for 
h = 1 and unit masses is 
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(i) 



where |A| < uo\UJ2 (attractive potential) and we allow for 
frequency diversity. The free Hamiltonian is diagonal- 
ized by a rotation in the x\ — X2 plane, with the angle 
that gives the eigenvectors a function of the cou- 

pling tan 20 = 2A/cc>2 — oj\. Master equations for both a 
common bath (CB) and separate baths (SB) have been 
compared by also analyzing entanglement decay time in 
Ref. [10], where both the similarity of the frequencies of 
the oscillators and the coupling strength were shown to 
contribute to the preservation of entanglement for a CB, 
leading to asymptotic entanglement in the case of identi- 
cal frequencies [3, [H, [l3| . The transition from SB to one 
CB underlies the capability of entanglement generation 
discussed in Ref. [14] and a physical implementation of 
the latter has been proposed recently [151 ]. 



Following Ref. [lOj, the system dynamics is described 
by a master equation that is valid in the weak-coupling 
limit between the system and environment, without the 
rotating- wave approximation [16]. Even if the obtained 
master equation has the same form as the exact one @, 
the coefficients are approximated for weak coupling 7 be- 
tween the system and environment. Taking this equation 
for strong coupling can lead to unphysical values for the 
reduced density and violation of positivity can appear at 
low temperatures and for certain initial states [lfj. In 
the following we restrict our analysis to small 7 = O.OIcji 
where we never encounter any unphysical dynamics. This 
is consistent with the fact that deviations of this mas- 
ter equation from one in the Lindblad form (preserving 
positivity) are in fact small for high temperatures (here 
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T = IOcji in natural units). Particularly useful for the 
purpose of understanding the physical behavior of the os- 
cillators dissipation is the master equation in the basis of 
the normal modes of the system Hamiltonian as given in 
Appendix [A] 

An important observation is that our results do not 
in fact depend on the specific choice of this master equa- 
tion. In particular, in Appendix[B]we compare our results 
with the one obtained from a master equation in the Lin- 
blad form, obtained by a rotating-wave approximation. 
Within this approximation the master equation is known 
to be in the Lindblad form [l6|, [I?} and we find almost 
exactly the same results as that obtained with the mas- 
ter equation (|A4|) . Therefore the phenomena predicted 
in the following do not depend on the specific details of 
the considered master equation. 



III. SYNCHRONIZATION 

The dynamical behavior of the two oscillators can be 
analyzed through their average positions, variances, and 
correlations, as we deal here with Gaussian states. The 
presence of a CB or of two (even if identical) SB leads 
to different friction terms in the dynamical equations of 
both first- and second-order moments [To[ with profound 
consequences. We recall that for a CB only the sum 
of positions x+ = x± + X2 is actually dissipating and 
this does not coincide with X + unless the oscillators are 
identical. 

Figure Ufa) shows the variance dynamics of two os- 
cillators starting from two vacuum squeezed states. To 
quantify synchronization between two functions f(t) and 
g(t), we adopted a commonly used indicator, namely, 



Cf, g (t,Ai) = Sf5g/ySf 2 Sg 2 where the bar stands for 

a time average / = J^ +At dt'f(t') with time window 
At and Sf = / — /. For similar evolutions \C\ ~ 1, 
while | C | ~ for different dynamics. The position 
variances for CB [Fig. QJa)] show a transient dynam- 
ics without any similarity between them, also in an- 
tiphase (C^Xa 2 ) < 0) 7 before reaching full synchroniza- 
tion [Fig.^b)]. 2 

A comprehensive analysis shows that this behavior is 
robust considering (i) different initial conditions, (ii) any 
second-order moments of the two oscillators (either of 
positions Xi^ or momenta pi^ or any arbitrary quadra- 
ture), and (iii) a whole range of couplings and detunings. 
Regarding (i), an important observation is that while in 
an isolated system the dynamics is strongly determined 
by the initial conditions, this is not the case in the pres- 
ence of an environment. After a transient (in which the 
initial conditions have an important role), we actually 
find synchronization independently of the initial state, 
with detuning and oscillator coupling therefore being the 
only relevant parameters. The full analysis (iii) for a CB 
allows us to conclude that synchronization arises faster 
for nearly resonant oscillators and that the deteriorat- 
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FIG. 1: (a) {x\{t))/S (light color) and (x 2 2 (t))/S (black) nor- 
malized with the shot noise <S, for 002 = \Au\ and A = 0.7lj 2 
starting from squeezed states for a CB and (b) synchroniza- 
tion C( x 2)( x 2) (£, 15) for a CB (dark color) and SB (light color) 
for temperature T = lOc^i (in natural units). The insets show 
synchronization values \C^ X 2^ X 2^\ varying UJ2/0J1 and X/uf at 
t = 300. Time t is scaled with uj\ (and is therefore adimen- 
sional) and 7 = O.OIcji, here and in the following. The initial 
state is separable with squeezing parameters 2 and 4, respec- 
tively, in the two oscillators. 



ing effect of detuning can be proportionally compensated 
for by strong coupling, as represented in the CB inset of 
Fig. lib). 

Moving now to the case of separate baths, a completely 
different scenario appears. The quality of the synchro- 
nization is generally poor (small |C|), not improving in 
time and dependent on the initial condition. The full pa- 
rameters map for \C\ is shown in the SB inset of Fig.QJb). 
In this case the oscillators do not synchronize in spite of 
their coupling even considering long times when finally 
the system thermalizes. 

The appearance of a synchronous dynamics only for 
a CB can be understood by considering the time evolu- 
tion of the second-order moments. The matrix govern- 
ing their time evolution [To| (see also Appendix [XJ) has 
complex eigenvalues {/^} (i = 1,..,10), named in the 
following dynamical eigenvalues. Their real and imagi- 
nary parts determine the decays and oscillatory dynam- 
ics of all second-order moments and variances. As shown 
in Fig. Efa), when A = all the eigenvalues are along 
the line —0.01 and for increasing coupling in the case 
of one CB they move in the complex plane assuming 
three different real values. In contrast, for SB all dy- 
namical eigenvalues have similar real parts that remain 
almost unchanged when varying parameters. Hence for 
SB the ratio between the maximum and the minimum 
eigenvalues Re(/iM)/Re(/x m ) ~ 1 is almost constant for 
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all parameters, while for a CB and for parameters for 
which synchronization is found [CB inset in Fig. OJb)] 
Re(/iM)/Re(/i m ) << 1 5 as shown in Fig. Efb). In this 
parameters regime, after a transient time, only the least- 
damped eigenvector survives, thus fixing the frequency of 
the whole dynamics of the moments. As a consequence 
of this mechanism, synchronization is observed by con- 
sidering the expectation values of any quadrature of the 
oscillators as well as higher-order moments. 



0.90_ 




1.00 



0.75 



FIG. 2: (a) Eigenvalues \ii in the complex plane for a CB for 
UJ2/0J1 — 1.31 and increasing the coupling from A = [circles 
with Re(/ij) rsj —0.01] to A = 0.9cji in the direction of darker 
colors, (b) Ratio between minimum and maximum eigenvalue 
Re(/x m )/Re(/iM) for a CB. 

We obtain an approximated analytical estimation of 
time scales by considering the master equation in the 
eigenbasis {X±,P±} of the free Hamiltonian [Eq. (pQ, as 
detailed in Appendix [XJ The master equations for com- 
mon and separate baths have the same expression in the 
case of detuned oscillators and the nature of dissipation 
(a CB or SB) appears only in the form of the damping 
coefficients. By eliminating the oscillating terms in the 
dynamics in the interaction picture, one obtains that, 
within this approximation, the decay rates of (P±) are 
given by 



f^SB 
1 ±± 



c 2 r n 



s 2 T 22 ± csT 12 



for SB and 

fCB _ 
1 ±± ~ 



(c ± s)(cr u ± s r 22 ) + (1 ± 2sc)ri 2 



(2) 



(3) 



for a CB, where c = cos 0, s = sin (with previously de- 
fined diagonalization angle of Hs), and rn } 22,i2 appear 
in the original master equation (see the Appendix of Ref. 
[Iq|). These approximated decays for the variances, to- 
gether with their average (f + f ++)/2 (for (P + P_)), 

for a CB and SB do agree very well with the real parts 
of the dynamical eigenvalues. 

As mentioned before, synchronization (for a CB) is 
observed by looking at the dynamics of both first- and 
second-order moments: the ratio between minimum and 
maximum dynamical eigenvalues is the same in both 
cases. Still, our interest is in the second-order moments 
due to their relevance in the quantum information shared 



by the oscillators. Furthermore, inspection of first-order 
moment dynamics allows us to establish connections with 
what is known in classical systems [1]. Two studied sce- 
narios for classical synchronization are the diffusive cou- 
pling where both oscillator dampings depend on the dif- 
ference of the velocity and the direct coupling where each 
one depends on the velocity of the other |l|. The quan- 
tum harmonic oscillators considered here for a CB display 
in their first-order moments a diffusive coupling up to a 
change of sign, which explains the antiphase character of 
their synchronization. 
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FIG. 3: (a) Mutual information and (b) discord on a log- 
arithmic scale for a common bath and separate baths (SB). 
Both the exact time evolutions and the filtered ones (Gaus- 
sian filter) are shown. The parameters are 0J2 — 1.05cji and 
A = 0.3cj?. 



IV. QUANTUM CORRELATIONS 

Once the conditions for synchronization to arise have 
been established, we explore this phenomenon by focus- 
ing on the information aspects, through mutual infor- 
mation shared by the oscillators and their quantum cor- 
relations. In particular, the total correlations between 
the oscillators are measured by the mutual information 
T(q) = S(qi) + S(q2) — S(g) where S stands for the von 
Neumann entropy, and Qi^ 2 ) is the reduced density matrix 
of each harmonic oscillator. A possible partition of corre- 
lations into quantum and classical parts that has received 
much attention recently is given by the quantum dis- 
cord S(q) = mm {Ui} [S(g 2 ) - S(g) + S{ 6l \{Ui})], 

with the conditional entropy defined as S(gi\{Uj}) = 
52iPiS(Qi\ni), Qi\Ui = lUo^h/Pi density matrix after a 
complete measurement {II j} on the second oscillator, and 
Pi = Tri2(n^£>). The importance of this measure of the 
quantumness of correlations relies on its capability to dis- 
tinguish and understand classical and quantum behaviors 
[20[ . Quantum discord has been reported recently also for 
continuous variables in Gaussian states [21] . 

Dissipation degrades all quantum and classical correla- 
tions [22|. Nevertheless, important differences are found 
when comparing a CB and SB for the same parameter 
choices. In Fig. [3] we show a fast decay of all the total 



[Fig. 13(a)] and quantum [Fig. [3(b)] correlations for SB. 
In contrast, for a CB we find that after a short transient 
both mutual information and discord oscillate around an 
almost constant value and their decay is nearly frozen. 
For these parameters and a common environment the os- 
cillators synchronize and = 0.95 at t ~ 270. The 
robustness of the quantum correlations for long times in 
synchronizing oscillators in a CB and the deep differences 
with the case of SB is surprising also because their respec- 
tive asymptotic values are really similar for detuned os- 
cillators. In other words, the upper CB curve in Fig. [3(a) 
[or Fig. [3(b)] will eventually thermalize, converging to a 
value very similar the one obtained for SB, while strong 
differences in the asymptotic values actually appear only 
in the case of identical oscillators 0- As a further result, 
the effect of increasing the temperature is mostly on the 
asymptotic state while the main features of the dynamics 
described here are still observed. 

We now focus on the case of a CB to look for spe- 
cific quantum features of the synchronization exploring 
different parameters regimes. The comparison of mu- 
tual information and discord in cases in which there is 
synchronization or the system dissipates without having 
time to synchronize is given in Fig. [4] (upper and lower 
curves, respectively) where we filter out the fast oscilla- 
tions to highlight the decay dynamics. For small cou- 
pling and large detuning, discord (shown in Fig. [4] for 
\/uo\ = 0.3 and 002/^1 = 1-4) and mutual information 
are rapidly degraded. In this case, when t = 200 there 
is no synchronous dynamics and C^ x 2^ x 2^ ~ 0. In con- 
trast, for strong coupling or for small detuning, synchro- 
nization occurs fast: for \juj\ = 0.8 and 002/^1 = 1-05, 
C( x \){xl){t — 200) ^ 1. In this case, after a short tran- 
sient, the dynamics of discord is almost frozen and it re- 
mains robust against decoherence. In exploring different 
parameter regimes we conclude that fast decay of classi- 
cal and quantum correlations is found in cases in which 
there is no synchronization while the emergence of syn- 
chronization accompanies robust correlations against dis- 
sipation (frozen decay). The inset in Fig. [4] represents the 
value of the discord after the fast decay (here for t = 300), 
where it is expected to be already in the plateau. There is 
a rather suggestive similarity to the synchronization map 
for a CB shown in the inset of Fig. [1(b). Considering that 
also entropy shows in this regime a slow growth, we con- 
clude that synchronized oscillators are characterized by 
a reduced leakage of information into the environment. 

One might wonder if the presence of a synchronous 
dynamics has any effect on entanglement, as in contrast 
with pure states mixed states with large quantum corre- 
lations can have even vanishing entanglement [23]. The 
presence of the environment for oscillators with different 
frequencies leads to a complete loss of entanglement in 
finite short times unless the couplings to the CB are "bal- 
anced" [Tcj. In the general case of detuned oscillators, 
even for large couplings, entanglement decay is typically 
faster than the time scales at which the system reaches 
synchronous dynamics both for a CB and SB, mostly at 




1000 



FIG. 4: Evolution of the discord for a CB and the parameters 
U2/001 = 1.05 and A = 0.8cji (dashed line A) and U2/001 = 
1.4 and A = 0.3cji (solid line B). The inset represents the 
quantum discord at t = 300 for a CB. 



this temperature (T = IOcji). Still, longer survival times 
for entanglement in a CB are found for small detunings 
and strong couplings [To[. 



V. DEPENDENCE ON INITIAL CONDITIONS 

We mentioned before that initial conditions do not play 
any important role in the appearance of synchronization. 
Indeed synchronous dynamics of the moments appears 
when an eigenmode dominates because of its slow dissi- 
pation rate and this goes beyond the specificity of the 
choice of the initial state. However, the details of the 
dynamics do depend on the latter as we illustrate for the 
following initial conditions: 

(i) the separable vacuum state: p = |0)(0| (g) |0)(0| 

(ii) the two-mode squeezed states 

P = tfi2(r)(|0><0|®|0)(0|)E4(r) 



where U\2(r) = exp 



and ai(a\) are 



i 2 — aia 2 )/2 

the usual annihilation (creation) operators; and (iii) the 
separable squeezed state: 

p = ^(r^lOXOl^Cn) (8) C/2(r 2 )|0><0|^(r 2 ) 5 

with Ui(ri) = exp -r(a\ 2 - af)/2 

Instead, quantum correlations S depend on the initial 
condition in the sense that more or less of the latter will 
be present. However after the short transient they always 
reach a plateau where information leakage to the environ- 
ment is greatly reduced. Both information leakage reduc- 
tion and synchronization are part of the same underlying 
phenomenon: that of a dissipation channel being much 
slower than the other. This behavior is seen in Fig. [5] 
where quantum correlations and the synchronization in- 
dicator are displayed for different initial conditions. We 
must further stress here that, since the asymptotic ther- 
mal state has S ~ 10 -4 , the plateau is expected to be 
very long. 
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FIG. 5: (a) Synchronization indicator and (b) decay of quan- 
tum correlations for different initial conditions in the case of 
a common bath: separable squeezed state with squeezings 
r\ — 2 and r<i — 4 (solid line), separable vacuum state (dotted 
line), and entangled two- mode squeezed state (dashed line) 
with squeezing r = 2. Other parameters are 002/(^1 = 1.1 and 
A = 0.8cj? 



CB and SB, the decay toward this value is clearly frozen 
in the presence of synchronization. In this case total and 
quantum correlations display a very slow decay (plateau) 
and the leak of information into the bath is reduced. 
The identification of the conditions for the occurrence 
of synchronization and its connection with quantum cor- 
relations reported here provides the path towards future 
extensions such as the study of arrays and networks, the 
analysis of the role of different environments, or the ex- 
ploration of eventual connections with biological systems 
in which synchronization is a widespread phenomenon. 
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VI. CONCLUSIONS 



Appendix A: Master equation in the normal modes 

basis 



Our analysis of the dynamics of dissipative quantum 
harmonic oscillators has allowed us to establish under 
which conditions synchronization appears. This phe- 
nomenon can appear in rather different forms but, in this 
paper it was reported during the transient dynamics of a 
(quantum or classical) system coupled to an environment 
and relaxing toward equilibrium. The emergence of syn- 
chronization was explained in terms of different temporal 
decays governing the system evolution and related to a 
separation between the eigenvalues of the matrix gener- 
ating the dynamics. We traced synchronization between 
second-order moment evolution from the existence of a 
slowly decaying eigenmode and found approximated ex- 
pressions for the variance decay coefficients in very good 
agreement with the real parts of the dynamical eigenval- 
ues. We found that synchronization arises in the pres- 
ence of a common bath, but not for separate ones, for 
any strength of the coupling between oscillators. It could 
be of interest to study the transition between such dif- 
ferent scenarios [l4|. The relevant parameters are X/uf 
and UJ2/0J1 since the dependence on initial conditions is 
actually weak. 

We have then characterized mutual synchronization 
from a quantum information perspective through mu- 
tual information and quantum discord exploring their dy- 
namics in different parameters regimes. We found that 
there is a signature of synchronization in the information 
shared by the oscillators. We have shown that discord 
and mutual information are more robust when the oscil- 
lators synchronize. In spite of the fact that after ther- 
malizing the asymptotic discord is negligible for both a 



The master equations describing the evolution of the 
reduced density matrix of the system of two different os- 
cillators, up to second order in the coupling strength, 
have been reported in Ref. [Toj| for both a common and 
separate baths. Here we provide these equations in the 
basis of the eigenmodes of the Hamiltonian. We stress 
that the exact master equation at all coupling orders has 
the same structure as ours, the difference being in the 
form of its coefficients 0, [l6[. For weak coupling this 
equation is a very good approximation to the exact one. 
Furthermore, in Appendix [B] we show that the full evolu- 
tion almost perfectly matches that of a master equation 
obtained by a rotating-wave approximation, the latter 
being always positive due to its Lindblad form. 



1. Common bath 

In the case of a common bath, we assume for the inter- 
action Hamiltonian between the system and the environ- 
ment the form Hf B = ^ k gk(bk + b\)x +1 where x+ = 
x\ + X2, bk(b\.) creates (annihilates) an excitation (with 
energy Q k ) over the kth. mode of the bath, and where the 
coupling coefficients g k are related to the density of states 
of the bath J(ft) through J(ft) = (gl/Sl k ) S(Q-Q k ). 
We assume an Ohmic environment with a Lorentz-Drude 
cut-off function, whose spectral density is 
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7T A 2 + f2 2 ' 
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The results represented in this work are obtained with 
bath-system coupling 7 = 0.01u;i and cutoff A = 50gji. 
The eigenvectors of H$ are 



X_ = cos 6x1 — sin #£2, 
X + = cos 6x2 + sin 9x\ , 



(A2) 
(A3) 



with tan 2(9 = 2A/cj| 



cies fl 
, ,2 1 , 



1 • 



The system eigenfrequen- 



are always different due to the coupling 2Vt 



uu 1 -r uJ% =t V^ 2 + (^2 — ^l) 2 - Neglecting energy renor- 
malization, the master equation in this eigenvector basis 
reads 



rln(t) D CB D CB 

-^1 = -i[H s ,p(t)] - [X_,p]] - -±±[X+, [X + ,p]] - D C + B [X + , [X-,p]] 

pCB pCB pCB pCB 

+-^[*_, [P-,p\] + -£M*+, [p+,p]} + -^[x+, [p-,p]] + -fM*-, [p+,p\] 
/ fcB f CB r CB r CB \ 

-i _^[X_,{P_,p}] + ^±[X + ,{P + ,p}] + ^[X + ,{P_,p}] + — ±[X_,{P + ,p}] , 



(A4) 



r 



with 



D\ 
D 

D 



CB 

CB 
++ 

CB 



(c - s)(cDn - sD 22 ) + (1 - 2sc)D 12 
(c + s)(cD n + sP> 22 ) + (1 + 2sc)D 12 



c 2 -s 2 



+- 



= 

T^CB _ 



1 ++ 
1 -+ 
1 +- 



(Ai + £22 + 2D12) + sc(Z>n - ^22) 



(c - s)(cF n - sF 22 ) + (1 - 2sc)F 12 
(c + s)(cPn + sP 22 ) + (1 + 25c)Pi 2 
(c - s)(cF 22 + aFn) + (c 2 - s 2 )Fi 2 
(c + s)(cP n - sP 22 ) + (c 2 - s 2 )F 12 
(c-s)(cr 11 -sT 22 ) + (l-2sc)r 12 

(c + s)(cr n + 5 r 22 ) + (1 + 2sc)r 12 

(c-s)(cT 22 + sr 11 ) + (c 2 -s 2 )T 12 



(A5) 



(c + s)(crn- s r 22 ) + (^-^)ri 2 

where c = cos 0, 8 = sin#, and the dissipation (I\j) 
and diffusion (DijJF\j) coefficients are defined in the 
Appendix of Ref.' [U specifically in Eqs. (A18)-(A20). 
The related set of equations of motion is 



d{XiXj) 1 



dt 
dt 



-({X^Pjj + iXj,^}) 



(A6) 



^({XuPjV + tfaXj,^})) 



2 

^CB 



(Vff + fffXPiPj) - miPjP^) 



CB 



^CB 



(PiP- 



D 



CB 



(A7) 



d({Xi,Pj}) 
dt 



)2/y.v.\ x.pCB 



- 3-3 \" ^ -J/ 1 ~M 

2(P i P i )-2^(X,X i ) , 
- tfJiiXuPjV-rfZjiXiP.j) (A8) 

where z,j = +, — . 

The time evolution of the vector R of all ten moments 
can be written in a compact matrix form 

R = M~R + N. (A9) 



The complex eigenvalues of M are {/^} with i = 1...10, 
elsewhere referred as dynamical eigenvalues. 



2. Separate baths 

When the two environments are thought to be iden- 
tical and independent from each other, the interaction 
Hamiltonian becomes 



(A10) 



1=1 k 



where the annihilation (creation) operators b % k (b % ^) be- 
long, respectively, to the zth thermal bath. The density 
of states of both of them is that of Eq. (|Aip . with the 
same 7 and A introduced before. The master equation 
has the same structure of Eq. (|A4ft but the coefficients 
are modified as follows: 



pSB 


= c 2 D n - 


f- s 2 D 22 - 


- 2csD 12 


pSB 


= c 2 D n - 


f- s 2 D 2 2 H 


- 2csD 12 


b SB _ 


= cs(D u 


- D 22 ) -{ 


-{c 2 -s 2 )D l2 


PSB 


= c 2 F n 4 


- s 2 F 22 - 


2csF 12 




= c 2 F„4 


- s 2 F 22 + 


2csF 12 


PSB 


= F™ = 


cs(F 11 - 


F 22 ) + (c 2 - s 2 


YSB 


= c 2 r n 4 


- s 2 r 22 - 


2csr i2 


i^SB 
1 ++ 


= c 2 r n 4 


- s 2 T 22 + 


2csT 12 


f^SB 
1 - + 


— 1 + _ — 


cs(T u - 


r 22 ) + (c 2 -s 2 ] 



Once these coefficients are used instead of those coming 
from a common bath, the equations of motion are for- 
mally identical to those of Eqs. dA6]) - (lA8]) . 
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Appendix B: Independent decay rates 

The eigenmodes (|A2|) and (|A3|) diagonalize the Hamil- 
tonian of the system Hs but are still indirectly cou- 
pled through the heat bath(s) as seen from Eqs. (|A6|) - 
(|A8p . This means that the eigenmodes cannot be consid- 
ered as independent channels for dissipation. Yet if we 
rewrite their master equation in interaction picture, we 
can neglect fast oscillating terms, as usual in the rotating- 
wave approximation, that is, eliminate exponents such as 
e ±i(n++n-)t ^ ue £ Q th e i r highly oscillatory behavior in 
comparison with the overall slower dynamics [13, HH • ^ 
we takethatwhich also rotate though more slowly, those 
such as e ±*(^+-^-) £ 5 an d keep only nonrotating terms. 
Finally, this procedure leads to an effective total decou- 
pling of the eigenmodes, which then dissipate indepen- 
dently to the heat bath(s) with the decay rates T±± and 
T±± [Eqs. (2) and (3)]. In some sense this time averag- 
ing approximation can be seen as renormalizing all dis- 
sipation coefficients having mixed indices +— (and — +) 
to zero, hence rendering the master equation as a ten- 
sor product of two independent evolutions. This could 
seem a bit too far fetched, but a comparison of the full 
dynamics and this approximation seems to be quite ac- 
curate as clear from Fig. [6j where we compare T±± and 
their average with the three different values of Re(/ii). 




-0.020 



0.0 0.1 



FIG. 6: Rates F 



CB 



from Eq. Q and and (f + 



T )/2 (dots) are compared with the (three different) real 



parts of the dynamical eigenvalues Re(/Xi 
in the case of a common bath for 002/(^1 = 



(continuous line) 
1.31. 



Inspection of these analytical expressions when varying 
system parameters confirms that dynamical eigenmode 
decays for SB all have similar real parts (in general F±2 
small and Tn ^ ^2)? while for a CB the decays can be 
significantly different (a factor 20 in Fig. [6]). This differ- 
ence between decay rates can be up to several orders of 
magnitude for parameters where synchronization appears 
faster. Synchronization is therefore linked to imbalanced 
dissipation rates of the eigenmodes, allowing the mode 
that survives longer to govern the dynamics. Within the 
discussed approximation its frequency is found to be 2^_ 



with Sl_ previously defined as the frequency of the eigen- 
mode X_ . This is independent of bath coefficients and 
we find very good agreement with the exact frequency. 

It can be easily seen that the rotating- wave approxima- 
tion described in this appendix, neglecting all highly os- 
cillatory terms with exponents e ± ' l ( Q + ±Q -) t ? leads to (CB 
and SB) master equations in the Kossakowski-Lindblad 
form [16]. In particular, in the case of a common bath it 
can be found that 



dp 
~dt 



-i[H s ,p(t)] (Bl) 
E + ^ \AAiP + pAAi - 2A iP A}\ 

E (DgBA y " IMP + P^A - 2A\p Al ] 



where A± = y / Q±/2X± + i/ \/ 2Q±P± . In spite of formal 
differences between Eqs. (jA4| and (jBl|) we actually find 
very good agreement between their dynamical evolutions. 
In Fig. Owe show that, in the limit of weak coupling con- 
sidered here, predictions for synchronization and discord 
are actually almost indistinguishable. Maximum devi- 
ations in this case are at least two order of magnitude 
smaller than the represented quantities. As expected, 
deviations increase for stronger system-environment cou- 
pling. The really weak dependence on the details of the 
master equation [Eqs. (|A4|) and (|B1|) ] when comparing 
the dynamical behavior of synchronization and quantum 
correlations between the oscillators further strengthens 
the generality of our results. 




FIG. 7: (a) Synchronization and (b) discord obtained from 
the complete master equation (|A4|) compared with the values 
obtained after the rotating-wave approximation as described 
in the text in the case of common bath for 002/^1 = 1.4, A = 



[1] A. Pikovsky, M. Rosenblum, J. Kurths, Synchronization: 
A Universal Concept in Nonlinear Sciences (Cambridge 



University Press, Cambridge, 2001). 



8 



[2] S.H. Strogatz, Nonlinear Dynamics And Chaos: With 
Applications To Physics, Biology, Chemistry, And Engi- 
neering, (Westview, BOulder, CO, 2001). 

[3] I. Goychuk, J. Casado-Pascual, M. Morillo, J. Lehmann, 
and P. Hanggi, Phys. Rev. Lett. 97, 210601 (2006); O. V. 
Zhirov and D. L. Shepelyansky, ibid. 100, 014101 (2008); 
S-B. Shim, M. Imboden and P. Mohanty, Science, 316, 
95 (2007). 

[4] M.Zhang, G. Wiederhecker, S. Manipatruni, A. Barnard, 

P. L. McEuen and M. Lipson. larXiv: 11 12.3636V 1. 
[5] G. Heinrich, M. Ludwig, J. Qian, B. Kubala and F. Mar- 

quardt, Phys. Rev. Lett. 107 043603 (2011). 
[6] C. A. Holm es, C. P. Meaney, G. J. Milburn, 

larXiv:1105.20861 
[7] J. P. Paz and A. J. Roncaglia, Phys. Rev. Lett. 100, 

220401 (2008); J. P. Paz and A. J. Roncaglia, Phys. Rev. 

A 79, 032102 (2009). 
[8] K-L. Liu and H-S. Goan, Phys. Rev. A 76 022312 (2007) 
[9] B. L. Hu, J. P. Paz and Y. Zhang, Phys. Rev. D 45 2843 

(1992); C-H. Chou, T. Yu, B. L. Hu, Phys. Rev. E 77 

011112 (2008) 

[10] F. Galve, G. L. Giorgi and R. Zambrini, Phys. Rev. A 

81, 062117 (2010). 
[11] T. Rocheleau, T. Ndukum, C. Macklin, J. B. Hertzberg, 

A. A. Clerk, and K. C. Schwab, Nature (London) 463, 

72 (2010); 

D. Van Thourhout and J. Roels, Nature Phot. 4 211, 
(2010); 

F. Marino, F. S. Cataliotti, A. Farsi, M. Siciliani de 
Cumis, and Francesco Marin, Phys. Rev. Lett. 104, 
073601 (2010); 

P. Verlot, A. Tavernarakis, T. Briant, P.-F. Cohadon, 
and A. Heidmann, Phys. Rev. Lett. 104, 133602 (2010); 
[12] K. R. Brown, C. Ospelkaus, Y. Colombe, A. C. Wilson, 



D. Leibfried, and D. J. Wineland, Nature (London), 471, 
196 (2011); 

M. Harlander, R. Lechner, M. Brownnutt, R. Blatt, and 
W. Hansel,^. 471, 200 (2011). 
[13] J. S. Prauzner-Bechcicki, J. Phys. A: Math. Gen., 37 
L173 (2004) 

[14] T. Zell, F. Queisser, and R. Klesse , Phys. Rev. Lett. 

102, 160501 (2009). 
[15] C. Cormick and J. P. Paz, Phys. Rev. A 81, 022306 

(2010) 

[16] H. P. Breuer and F. Petruccione The Theory of Open 
Quantum Systems (2002, Oxford University Press). 

[17] A. Rivas, A. D. K. Plato, S. F. Huelga and M. B. Plenio, 
New J. Phys. 12, 113032 (2010) 

[18] H. Ollivier and W. H. Zurek, Phys. Rev. Lett. 88, 017901 
(2001). 

[19] L. Henderson and V. Vedral, J. Phys. A 34, 6899 (2001). 
[20] A. Datta, A. Shaji, and C. M. Caves, Phys. Rev. Lett. 

100, 050502 (2008); B. P. Lanyon, M. Barbieri, M. P. 

Almeida and A. G. White, ibid. 101, 200501 (2008); L. 

Mazzola, J. Piilo, and S. Maniscalco, ibid. 104, 200401 

(2010); T. Werlang, C. Trippe, G. A. P. Ribeiro and G. 

Rigolin, ibid. 105, 095702 (2010). 
[21] P. Giorda and M. G. A. Paris, Phys. Rev. Lett. 105, 

020503 (2010); G. Adesso and A. Datta, ibid. 105, 030501 

(2010). 

[22] G. L. Giorgi, F. Galve, and R. Zambrini, Int. J. Quantum 

Inf. 9, 1825 (2011). 
[23] F. Galve, G. L. Giorgi, and R. Zambrini, Phys. Rev. A 

83, 012102 (2011). 
[24] R. Vasile, S. Olivares, M. G. A. Paris, and S. Maniscalco, 

Phys. Rev. A 80 062324 (2009). 



